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The vortex lattice (VL) in the high-re superconductor YBa2Cu3C>7, at 2 K and with the magnetic 
field parallel to the crystal c-axis, undergoes a sequence of transitions between different structures 
as a function of applied magnetic field. However, from structural studies alone, it is not possible to 
determine precisely the system anisotropy that governs the transitions between different structures. 
To address this question, here we report new small-angle neutron scattering measurements of both 
the VL structure at higher temperatures, and the field- and temperature-dependence of the VL form 
factor. Our measurements demonstrate how the influence of anisotropy on the VL, which in the- 
ory can be parameterized as nonlocal corrections, becomes progressively important with increasing 
magnetic field, and suppressed by increasing the temperature towards T c . The data indicate that 
nonlocality due to different anisotropies play important roles in determining the VL properties. 



PACS numbers: 74.25. Uv, 74.72. Gh, 61.05.fg 
I. INTRODUCTION 

In an applied magnetic field, type-II superconductors 
can be investigated directly via studies of the magnetic 
vortex lattice (VL). One of the most direct experimen- 
tal probes of the VL is small-angle neutron scattering 
(SANS). Being a bulk sensitive technique, SANS con- 
tinues to play a pivotal role in clarifying the details 
of the mixed state in various classes of superconduc- 
tors^— Typically, in SANS studies it is reported how 
the VL structure and coordination vary as functions of 
both magnetic field and temperature (T). The results 
are often used to provide evidence as to the sources of 
anisotropy that might cause the VL structure to deviate 
from the perfectly hexagonal coordination predicted for 
an isotropic material^ By increasing the applied mag- 
netic field, these anisotropies become increasingly influ- 
ential on the inter- vortex interaction, and can cause tran- 
sitions between different VL structures. 

The simplest theoretical models that can be used to 
predict the occurrence of VL structure transitions are 
those based on local electrodynamic theory ('London' 
theory) extended with nonlocal correction terms^ These 
correction terms provide a direct coupling between a 
system anisotropy and the VL properties. The size of 
the correction terms scale according to re~ 2 , where re is 
the Ginzburg-Landau parameter. Hence for high-re ma- 
terials such as cuprates, nonlocal corrections are often 
neglected and local theory is deemed sufficient to un- 
derstand observation. However, VL studies of reason- 



ably high-re materials, such as borocarbides^ and s-wave 
VaSij^ show that in these materials it is necessary to 
consider nonlocal effects associated with a Fermi surface 
anisotropy in order to understand the field-dependent VL 
structure^-— Furthermore, in high-re materials that ex- 
hibit a d-wave gap anisotropy, nonlocal corrections can- 
not be neglected because the nonlocal length-scale £o is 
momentum-dependent £o(k) oc 1/Ak and hence will be 
divergent at the nodes. As a direct consequence of these 
'<i-wave' nonlocal effects, unconventional VL structures 
are predicted*^— More sophisticated models can also 
account for the effects due to anisotropy suggested by the 
aforementioned nonlocal theories. Within quasiclassical 
Eilcnberger theory, both the Fermi surface and super- 
conducting gap anisotropies can be considered in equal 
measure when evaluating the VL free energy, and again 
unconventional VL structures are expected as functions 
of both field and Ti 17 i 21 ~ — Therefore, while nonlocal ef- 
fects are often neglected in high-re materials, there is sub- 
stantial evidence to suggest that they play a decisive role 
in determining the VL properties. 

Our subject high-re superconductor is YBa2Cu307_,5, 
the mixed state of which has been under scrutiny more 
than twenty yearsi^ir— SANS studies in particular have 
made important contributions towards understanding 
the mixed state in this material^-— with the motiva- 
tion to understand the physics of both the VL and the 
pairing state within an anisotropic host material that has 
a nodal gap symmetry. YBa2Cu307_a exhibits an or- 
thorhombic crystal symmetry and is composed of stacked 
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two-dimensional Cu02 plane layers, and one-dimensional 
CuO chain layers. For optimally- and over-doped com- 
positions, the CuO chains display both long-range or- 
der along the crystal b-axis, and metallic behavior4i As 
a consequence it is expected that, in addition to Cu02 
plane states, electronic states primarily associated with 
the chains are also superconducting below T c . This pic- 
ture is supported by reported values of the in-plane pen- 
etration depth ratio, 7a (= A a /Ab oc \/m*/ml) lying in 
the range 1.2-1.5i 29 ' 37 i 40 ' 42 ~ — While the physical origin 
of the chain pairing interaction remains to be definitively 
explained, current evidence suggests that it may origi- 
nate from a proximity effect coupling between chain and 
plane states^ 5 - - — A further consequence of the crystal 
orthorhombicity is that the predominantly d x 2_ y 2 order- 
parameter must contain an additional and finite s-wave 
admixture? 4 ^ Evidence for this admixture is provided by 
phase sensitive, 49 tunneling, 50 and /uSR 5 ^ studies. 

The striking result of earlier SANS studies on fully- 
doped, albeit crystallographically twinned YBa2Cu307, 
was the observation that with increasing magnetic field, 
the VL structure evolves continuously from a low field 
hexagonal structure towards an almost square structure 
by 11 T i 37 ' 39 The physical origin of this transition was 
suggested to be due to the increased importance of the 
d-wave order parameter anisotropy as vortex cores move 
closer together at high field i 17 ' 22 ' 23 However, it was also 
acknowledged that such a transition may also be driven 
by nonlocal effects associated with a band structure 
anisotropy. 37 ' 39 More recently, we reported a SANS study 
of the VL in a sample of twin-free, and fully oxygenated 
YBa2Cu30y4S In this sample the pinning effects on the 
VL due to twin boundaries, a characteristic feature of 
all previous SANS studies^ - — are almost entirely sup- 
pressed. For the first time, this allowed a clear observa- 
tion of the intrinsic VL structure, and we reported an 
entirely new sequence of field-dependent VL structures 
and structure transitions at 2 K, with all of the structure 
phases separated by first-order transitions^ - 

In this paper, we extend our previous SANS study of 
twin-free YBa2Cu307 with new measurements of both 
the T-dependent VL structure, and the field- and In- 
dependence of the VL form factor. Measurements of 
the form factor have proved to be remarkably useful 
in characterizing the mixed state of other superconduc- 
torw£— — — and are often used in conjunction with sim- 
ple theories in order to extract measures of the penetra- 
tion depth and coherence length. The results reported 
here represent the first detailed investigation of the VL 
form factor in YBa2Cu307. 



II. EXPERIMENTAL METHOD 

The SANS investigations reported in this study have 
been carried out using the SANS-I and SANS-II instru- 
ments at the Paul Scherrer Institut, Villigen, Switzer- 
land, the Dll instrument at Institut Laue-Langevin, 



Grenoble, France, and the SANS-NG3 instrument at the 
NIST Center for Neutron Research, Gaithersburg, USA. 
For all experiments, neutrons of a mean wavelength be- 
tween 6 and 10 A were selected with a 10% FWHM 
spread in wavelength. The neutrons were collimated over 
distances between 6-14 m before the sample position, 
and the diffracted neutrons were collected by a position- 
adjustable two-dimensional multidetector. 

All of the SANS measurements were performed using 
the same twin-free sample of YBa2Cu307 that was also 
used in Ref.fi^. A mosaic of six single crystals was assem- 
bled with a total mass of ^30 mg. Each single crystal was 
grown from flux within home-made BaZr03 crucibles. 56 
The as-grown crystals were individually detwinned un- 
der uniaxial stress at elevated T ; 57 ' 58 and subsequently 
oxygenated under a high pressure oxygen atmosphere 
to achieve the overdoped O7 phased The zero-field T c 
of the largest of the six crystals was measured using a 
SQUID magnetometer. The mid-point of the transition 
was found to be 89.0(5) K which is consistent with the 
sample being maximally dopedi^ - The sample mosaic was 
assembled on a 1 mm thick Al plate with the crystal ef- 
faces flat and co- aligned about the c-axis. X-ray Laue 
measurements showed the mosaicity of the co-alignment 
to be <1.5°. The detwinning ratio of the mounted sam- 
ple was determined from single crystal neutron diffrac- 
tion measurements of the (100) and (010) nuclear peak 
intensities. No signal attributable to a minority domain 
was detectable within the sensitivity of the measurement, 
limiting the minority domain population to < 1 %. 

The sample mosaic was loaded into the variable- 
temperature-insert of a horizontal field cryomagnet with 
the a-axis vertical, and such that the c-axis was paral- 
lel to the field direction and also approximately parallel 
to the neutron beam. In our previous report^ - we de- 
scribed two different techniques for the preparation of 
the VL in the sample. The first is the standard field- 
cool (FC) procedure. In this case, the VL is prepared 
by applying the field above T c , and then field-cooling 
through T c to the measurement T. The second is the 
oscillation field-cooling (OFC) procedure. In this case, 
the field-cool through T c is carried out in a field oscil- 
lating around the target value B, with an amplitude of 
typically 0.1 % to 0.2 % of B, and with an oscillation 
frequency of « 2 ruin -1 . We identified^ - the role of the 
OFC procedure as that which equilibrates the VL struc- 
ture into the preferred structural free energy minimum 
T, that exists for an effective T which is lower than the 
irreversibility temperature (T; rr ) exhibited when using a 
standard field-cool. The OFC procedure is observed to 
be most effective at low fields, where the stabilization of 
the intrinsic VL structure is further accompanied by sig- 
nificant improvements in VL quality. For these reasons, 
in this paper we only report SANS measurements of the 
optimized VL prepared by the OFC procedure. 

When at the desired measurement T, the OFC proce- 
dure is halted and the field is held constant at the target 
value. The subsequent diffraction measurements are car- 



3 




FIG. 1. (Color online) VL diffraction patterns obtained at 2 K, and in fields applied parallel to the crystal c-axis, of (a) 0.5 T, 
(b) 5.0 T and (c) 7.5 T. The axes indicated in (a) apply to all figures. In each image, dashed line patterns indicate the VL 
structures. Solid lines represent the basis vectors of the primitive cell, with the primitive cell opening angle of (a) <f>, (b) p 
and (c) v. The images are constructed by summing together the diffraction patterns obtained at a series of angles about the 
horizontal and vertical axes. This allows the presentation of all the Bragg spots in a single picture. Statistical noise at the center 
of the patterns has been masked, and the data smoothed with a Gaussian envelope smaller than the instrument resolution. 
The VL in real space can be visualized by rotating the reciprocal VL by 90° and adding an additional vortex at the origin. 



ried out by rotating the sample and cryomagnet together 
about the horizontal and vertical axes to angles that sat- 
isfy the Bragg condition for a certain Bragg spot on the 
detector. In all cases, background measurements were 
taken at T > T c , and subtracted from low T foregrounds 
to leave just the VL signal. Data visualization and anal- 
ysis was performed using the GRASP analysis package. 61 



III. RESULTS AND DISCUSSION 
A. Vortex Lattice Structure at 2 K 

1. Field-Dependence at 2 K 

To set the stage for the new results in this paper, we 
begin by reviewing our previous study.— At 2 K, and 
within the field range up to 10.8 T, we observe a se- 
quence of first-order VL structure transitions between 
the two orientations of the VL primitive cell permitted by 
the orthorhombic crystal symmetry. 62 Both of these ori- 
entations as shown in the new SANS diffraction patterns 
presented in Fig.[TJ Fig.[T](a) shows a diffraction pattern 
obtained in an applied field of 0.5 T. The overlaid hexag- 
onal pattern indicates that the VL structure is composed 
of a single distorted hexagonal domain aligned with the 
underlying atomic lattice. We refer to this structure type 
as the low field structure (LFS). The shape of the VL unit 
cell may be described by the angle <fi which lies between 
the two basis vectors indicated in Fig. [T] (a) . On increas- 
ing the field, the VL structure undergoes a first-order 90° 
re-orientation transition into another single domain dis- 
torted hexagonal structure. Fig. Q] (b) shows an example 
of this structure type in a diffraction pattern obtained 
at 5.0 T. We refer to this structure type as the interme- 
diate field structure (IFS). The sign of the distortion of 
the IFS remains consistent with that of the LFS; both 



structures are stretched from an isotropic hexagonal co- 
ordination along the a* direction. The overlaid pattern 
in Fig. [T] (b) indicates how we describe the primitive cell 
of the VL with a different characteristic angle p. On fur- 
ther increase of the field, we observe another first-order 
transition in the VL structure, this time between the IFS 
and a high field structure (HFS) that is dominated by 
a rhombic shape composed of four intense Bragg spots. 
This rhombic shape is highlighted in Fig. [1] (c) which 
shows an example diffraction pattern obtained from the 
HFS in an applied field of 7.5 T. Since the primitive cell 
of the HFS is of the same orientation as that for the LFS, 
in order to distinguish between them we use the symbol 
v to label the opening angle of the HFS primitive cell. 

An interesting detail of the diffraction pattern shown in 
Fig. [2(b), and characteristic of the IFS phase, is that the 
horizontal spots with q parallel to the b* axis (q || b*) 
are weaker in intensity than the other four spots. This 
is converse to the general expectation that weaker spots 
appear at larger q than more intense spots and are conse- 
quently considered to be higher-order spots of the recip- 
rocal VL. In the HFS phase, the more usual situation is 
realized; the weaker spots with q parallel to the a* axis 
(q || a*) exhibit a larger q than the four strong spots 
that make up the rhombic structure, and thus may be 
considered as higher-order spots. According to this de- 
scription, we see that the rhombic shape is distorted away 
from an isotropic square lattice by a stretching along the 
b* direction. Such a distortion is of opposite sign to that 
observed for the hexagonal structures; a possible physical 
reason for this is discussed later. Alternatively, the HFS 
may also be considered as a distorted hexagonal structure 
of the same orientation as the LFS phase, and the same 
direction of distortion as the LFS and IFS phases. While 
there is no formal method to distinguish between such 
descriptions of the VL structure, other data presented 
later will show that with increasing field, the spots of the 
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FIG. 2. (Color online). The field-dependence at 2 K of the 
primitive cell opening angle for the various VL structures. 
Inset schematics indicate the reciprocal VL primitive cell for 
each phase. The dashed line at 90° indicates the opening 
angle of a perfectly square VL. The dashed-dot line at 60° 
indicates the opening angle of an isotropic hexagonal VL. The 
(dark) gray shaded portions indicate the field regions over 
which we observe first-order VL structure transitions. The 
error bars are sufficiently small to be masked by the data 
symbols. Part of these data were first presented in Ref. l40l . 



rhombic structure become further dominant in intensity 
over the weak spots with q || a*. 

In Fig. [2J we show a quantitative analysis of the VL 
structure and plot the field-dependence of the primitive 
cell opening angle at 2 K. At 2.1(4) T and 6.7(2) T we ob- 
serve clear discontinuities in the overall field-dependence, 
confirming the first-order nature of the VL transitions. 
The widths in field of the regions of structure co-existence 
are estimated to be 0.75(5) T for the low field transition, 
and 0.45(5) T for the high field transition. The finite 
width of these co-existence regions may arise from slight 
differences between the various crystals of the sample mo- 
saic or from hysteresis within them, though the broader 
LFS - IFS transition also indicates that the LFS and IFS 
types are closer in free energy than the IFS and HFS 
phases. Fig. [2] also shows that the primitive cell angle of 
the HFS is still evolving with field until the highest fields 
available. Whilst at no field do we observe a perfectly 
square coordination with a primitive cell angle of 90°, 
we stress that in orthorhombic YBa2Cu30y there is no 
particular symmetry reason for the VL to stabilize into 
a perfectly square structure at high fields. 

Next we examine the distortion of the VL structure in 
the LFS and IFS phases. The stretching of the diffraction 
patterns along a* is consistent with that expected due to 
chain enhanced superconductivity along the b direction. 
This anisotropy is incorporated into local London theory 
by inserting a second-rank effective mass tensor into the 
expression for the London free energy^ - — For the uni- 
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FIG. 3. (Color online). The field-dependence at 2 K of the pa- 
rameter r) determined for both the LFS and IFS. The (dark) 
gray shaded portion indicates the field region over which we 
observe a first-order transition between the LFS and IFS 
phases. The line corresponds to a linear fit of the data ob- 
tained in the IFS phase. Error bars not visible can be consid- 
ered of order the size of the symbol. Part of these data were 
first presented in Ref. H^. 



axial case where m* = m£, 7a = \J™* a jm\ = 1) the VL 
structure is undistorted and perfectly hexagonal. In such 
a case, the spot distribution can be overlaid by a circle. 
For the bi-axial case m* ^ m£, the VL structure in re- 
ciprocal space is scaled by 7^ 2 along a* and by 77 
along b* , resulting in an elliptical spot distribution^^ 
In the anisotropic London limit, the axial ratio of the 
ellipse, which we call 77, gives directly the value of the 
anisotropy parameter 7a. Geometry shows that for the 
LFS, 77 can be determined directly from the defined open- 
ing angle <j> according to 7]lfs = \/3/tan(90° — {4>/2)). 
For the IFS, 77 can be found from the angle p accord- 
ing to tjifs = (\/3tan (p/2)) _1 . Fig. [3] show the field- 
dependence at 2 K of 77 obtained within both the LFS and 
IFS phases. For the LFS phase, the value of 77 remains es- 
sentially field-independent. The mean value taken across 
this phase is r\ = 1.28(1), a value similar to those found in 
previous studieS) 29 ' 37 ' 42 i 44 some of which were performed 
on YBa2Cu307_5 samples of slightly different composi- 
tion. On moving between the LFS and IFS phases, the 
first-order change in the VL structure is accompanied by 
a discontinuity in 77. Within local London theory, this 
would imply a similarly sharp discontinuity in the value 
of 7a which is not expected. This suggests that nonlocal 
effects are becoming important. On further increase of 
the field, the value of 77 is observed to fall monotonically, 
which similarly within the local theory would indicate a 
suppression of the in-plane anisotropy. However, results 
in the following sections will show that such a simple 
picture may be inappropriate, and that a more detailed 
theory is required in order to explain the observations. 
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2. Dtscussion 

Within the literature, a broad range of theoretical work 
attempts to describe the field-dependent behavior of the 
low T VL structure in high-K materials. In our case, the 
simplest model is the anisotropic London theory. This 
is a purely local theory that is valid for k» 1 and for 
fields B <C B C 2, a regime deemed suitable for experiments 
on High-T c materials J£r£2. According to local theory, the 
preferred VL structure is expected to be hexagonal in co- 
ordination, but with no preferred structural orientation. 
Since this orientational degeneracy is expected within lo- 
cal theory, even an extremely weak interaction may lead 
to a preferred orientation. Therefore local theory may 
still be dominant at low fields, even though preferred 
orientations are seen. Since no field-driven VL structure 
transitions are expected from the predictions of local the- 
ory, this demonstrates that the local approach neglects 
important details, even if the theory is able to describe 
the VL distortion in high-K materials. 

An orientational T can be found at every field and T 
after extending the local London theory with nonlocal 
corrections that take the form of higher-order terms in 
the expression for the London free energy. These higher- 
order terms provide a coupling of the VL both to the 
anisotropy of the band structure ^ 17 ! 19 and, in general, 
an anisotropy of the gapJ^~— In each case, numerical 
evaluations of the free energy are able to reproduce field- 
driven and first-order re-orientation transitions between 
hexagonal VL structures that are analogous to our ob- 
served LFS - IFS transition. However, nearly all model 
predictions are made using an idealized fourfold symmet- 
ric system, and predictions for twofold symmetric sys- 
tems are comparatively fewi 15 i 16 ' 66 Seemingly the most 
suitable predictions for twofold symmetric systems such 
as YBa2Cii307 are provided by the nonlocal theory of 
Kogan et al. 15 ' 16 However, before we discuss the applica- 
tion of the theory to our resultS ) 15 i 16 we next discuss the 
theory in more detail, because we have discovered it to 
possess some previously unreported properties. 

Kogan's nonlocal theory^ accounts for the effects of 
Fermi-surface anisotropy, but not gap anisotropy, on the 
VL structure (although we expect that a <i-wave gap 
anisotropy will have similar effects^ 9 -). The theory is for- 
mulated by developing the lowest-order nonlocal correc- 
tions to the local London relationship between supercur- 
rcnt density, j and magnetic vector potential, a. Kogan et 
al. showed that by including these corrections the general 
form of the nonlocal London equation becomes^— 
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where, rriij is a normalized effective mass tensor and 
A is the geometrically averaged penetration depth = 
(AajA^Az)" 1 / 3 . Nonlocality is captured by the term de- 
pendent on both the reciprocal VL vector G, and the 
fourth tank tensor, riiji m . This tensor provides a cou- 
pling between the VL supercurrents and the microscopic 
properties of the Fermi surface, and is dependent on 



the fourth moment of the Fermi velocity according to 
niji m oc (viVjViVm) J^r— Without the term dependent on 
riiji m , Eq. [T]is the local anisotropic London equation. 

The free energy density, F of a certain VL coordination 
is computed numerically by summing the contributions 
due to the various Fourier components of the field distri- 
bution described by Eq. [1] At each field, the preferred 
VL coordination is that which minimizes F. For fields 
applied along a fourfold symmetric z-axis, the x and y 
directions are equivalent, and F is evaluated by perform- 
ing the sum>i£ 
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where n\ oc (v x ) = (v y ) and n 2 oc (v 2 v 2 ). When evalu- 
ating F, the logarithmic infinity in the magnetic field at 
a vortex core given by London theory has to be regular- 
ized. This is done by inserting a 'core-correction' term of 
the form exp (— G 2 £ 2 ) into the numerator of Eq. [2j and 
may itself be anisotropic due to effective mass anisotropy. 
Since effective mass is represented by a second rank ten- 
sor, for the fourfold case, the effective mass of the carri- 
ers is isotropic in the basal plane and the effects of the 
core correction al 16 are thus taken to be isotropic too. 
Similarly the local limit for the value of the penetration 
depth is also isotropic since it too depends on the effec- 
tive mass. Therefore, the nonlocal effects are the only 
cause of anisotropy. 

For fields applied along a twofold symmetric z-axis, 
(appropriate for the case of fioH || c in YBa 2 Cu307) 
the effective mass is anisotropic (the x and y directions 
are inequivalent) so both the value of penetration depth 
and the core correction are modified to take this into ac- 
count^ Kogan et al. showed that the general expression 
for F becomes: 15 
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where to, is the effective mass in the zth direction, and 
the tensors n% oc (v x ), n\ oc (v y ) and 714 oc (v 2 v 2 ). 

By using Eqs. and 02 Kogan et al. calculated the 
field evolution of the VL structures for both the four- 
fold and twofold symmetric casesJ^ Their reported pre- 
dictions are shown in Fig. |4j For the fourfold case, the 
hexagonal VLs predicted in the low- and mid-field ranges 
break the fourfold crystal symmetry. Hence at any field 
where hexagonal VLs are expected, there are two de- 
generate VLs rotated 90° to one other. The lowest-field 
hexagonal VLs are predicted to be oriented at 45° to 
the mid-field hexagonal VLs, and change from one to the 
other via a first-order transition. The mid-field VLs then 
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FIG. 4. The predicted field-dependences of the VL opening 
angle f3 for both the fourfold symmetric case (down pointing 
triangles and diamonds), and the twofold s ym metric case (up 
pointing triangles). Figure taken from Ref Ha. 



distort with increasing field before locking via a second- 
order transition into a single square VL of definite orien- 
tation relative to the crystal axes. 

For the twofold symmetric case no low-field hexagonal 
VL expected, and the analogue of the mid-field hexagonal 
VL is present in only one orientation, aligned with one of 
the crystal axes. This VL is predicted to distort as the 
field is increased and to lock-in via a second-order transi- 
tion to a non-square high-field VL with field-independent 
shape, and again in a definite orientation relative to the 
crystal axes. Whereas the fourfold symmetric version of 
the theory gives reasonable agreement with observations 
in superconductors with appropriate field directions 
it does not agree with those obtained with the field along 
a twofold symmetric axis^i 

We have therefore reconsidered the twofold symmet- 
ric version of the theory, and discovered that it has 
some previously unreported properties. Starting in the 
twofold symmetric case, we can make a scale transfor- 
mation of the x- and y-components of the spatial vari- 
ation of the magnetic field: m x (G x ) 2 —> (G' x ) 2 and 
m y (Gy) 2 — > (G' y ) 2 , and thus remove the effective mass 
anisotropy from the theory. Now all the local and nonlo- 
cal contributions to the penetration depth become iden- 
tical in form to the fourfold symmetric case (with coeffi- 
cients derived from the twofold parameters also rescaled). 
This k-space scaling is also exactly that required to make 
the core corrections in the theory isotropic. Hence, the 
twofold symmetric version of the theory has exactly the 
same structure as the fourfold symmetric version, except 
that the scale transformation distorts the shape of the 
predicted VLs and re-scales the value of the field. 

Using Eq.[31 we have made a numerical implementation 
of the theory for the twofold case in order to study the 
consequences of the isomorphism of the twofold and four- 
fold theories. Firstly, we find that in the mid-field region, 
two degenerate hexagonal VLs are expected [Fig. [5], in- 
stead of just the single hexagonal VL (at low /3) that was 
identified previously. 15 This degeneracy may be expected 
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FIG. 5. Numerical evaluations of the free energy carried out 
using the Kogan nonlocal theory in the mid-field region and 
for the twofold symmetric case-^ Inset diagrams show the 
degenerate reciprocal VLs. Filled symbols denote the prim- 
itive cell, and empty symbols represent the additional spots 
required to complete the hexagonal structure. Dashed lines 
indicate horizontal and vertical symmetry planes. The param- 
eters used for the calculations are the same those in Ref. Ha. 
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FIG. 6. The field-dependence of the primitive cell opening 
angle /3 as calculated using the Kogan theory for the twofold 
symmetric case. In the mid-field region, we plot only the 
small-/3 solution in order to compare with Fig. [3] The recip- 
rocal VL structures are defined by the inset diagrams. The 
short- and long-dashed lines indicate the crystal symmetry 
planes for the fourfold case. For the twofold case, only the 
short-dashed line symmetry planes remain. The parameters 
used for the calculations are the same as those in Ref. [La. 



for the twofold case since for the fourfold symmetric case, 
two identical (but rotated by 90°) distorted hexagonal 
VLs have the same free energy. Secondly, we find that 
the minimum free energy at the lowest fields corresponds 
to VLs in orientations that are rotated approximately 45° 
from the mid-field VLs as shown in Fig[SJ^ After scaling 
to go from the fourfold to the twofold symmetric case, 
these lowest field VLs are in unsymmetrical orientations 
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relative to the crystal axes, and so were not searched for 
in the earlier work. Again, in agreement with the four- 
fold symmetric case, the low- and mid-held VL structures 
are separated by a hrst-order reorientation transition of 
the primitive cell. In the high-field region, the expected 
single non-square VL solution corresponds to a fourfold 
square VL distorted by the scaling factors. 

We also observe that there is no change of symmetry 
expected at the transition between the mid-field hexag- 
onal VL and the high- field non-square VL [Figs. HJ [6]. 
This is of concern since in a real twofold symmetric sys- 
tem there cannot be a second-order transition between 
such VLs: it is only predicted in the twofold case be- 
cause of its 'hidden fourfold symmetry'. This symmetry 
is a property of the theory and not a necessary property 
of the system being described. Hence, in practice, any 
property of the system not contained within the theory 
will remove the sharp transition in the twofold case, and 
replace it with a smooth crossover. We can remove this 
hidden symmetry easily by putting a different effective 
mass anisotropy into the core correction from that in the 
penetration depth. This maintains the twofold crystal 
symmetry, but removes the isomorphism of the two- and 
a fourfold theories. Once this is done, we find that both 
the the high-held second-order transition is removed, and 
that the degeneracy between the two mid-field hexagonal 
VLs is removed so that only one orientation is expected. 

We see that the predictions of the Kogan nonlocal the- 
ory depend sensitively on the relationship between the 
anisotropy of the cores and that of the penetration depth. 
In YBa2Cu307, we may expect the coherence length of 
the carriers localized in the CuO chains to be quite differ- 
ent in value and degree of anisotropy from the properties 
of the carriers in the Cu02 planes. In addition, Kogan's 
theory contains no gap anisotropy, and the VL core struc- 
ture in YBa2Cu307 must be affected by the positions 
of the gap nodes, which will also introduce a nonlocal 
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We can be fairly sure that these effects 



will remove the degeneracy between the two different ori- 
entations of distorted hexagonal VL that are predicted 
at each field for the twofold symmetric case. Our experi- 
mental results show that these two different orientations 
exist at fields above and below ^2.1 T [Fig[2]; it is possi- 
ble that a field-dependence of the core anisotropy could 
lead to the observed field-dependent transition between 
them. Nevertheless, in order to apply the twofold sym- 
metric form of Kogan's non-local theory to observations 
real materials, further extensions, perhaps in the form of 
additional higher-order anisotropy terms, are required. 
This conclusion is strengthened by the fact that we do 
not observe the very low field ^45° rotated VL (which is 
apparently expected for all values of anisotropy 17 ). Also 
the transition to our high-field phase is first-order, again 
in disagreement with simple nonlocal theory^ 

Despite the noted shortcomings of the Kogan theory 
in its current form, we nevertheless emphasize certain 
details within the predictions that may be compatible 
with our data. Firstly, we have revealed that both a 
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FIG. 7. (Color online) The zero field in-plane anisotropies 
of superconducting YBa2Cu3C>7 that are expected to influ- 
ence the VL properties, a) A schematic of the irreducible 
quadrant of the first Brillouin zone showing the main Fermi 
surfaces at k z = 0.— i 7 — b) The predominantly d x i_ y i order 
parameter combined with the s-wave admixture suggested by 
phase-sensitive measurements.— 



low field and first-order VL reorientation transition may 
be expected for the twofold symmetric case. After ap- 
propriate modification of the theory, it may be turn out 
that this transition will involve a 90°, and not a ~45°, 
re-orientation of the primitive cell, and thus be more 
directly analogous to our observed transition between 
the LFS and IFS phases. Secondly, we find that the 
high field VL is expected to evolve towards a distorted 
square-like structure that is analogous to the behavior 
observed in the HFS phase. While we have clarified that 
in a real twofold symmetric system the predicted high 
field second-order transition can not occur, the impor- 
tant detail is that at high field the vortex nearest neigh- 
bors are expected to approach directions parallel to min- 
ima in the Fermi velocity^ This observation is compat- 
ible with calculations of the zero field Fermi surface of 
YBa 2 Cu 3 07^iii shown schematically in Fig. (a). The 
dominant bands are large square-like contours associated 
with the Cu02 planes that are distorted so that the Fermi 
velocity is larger along k x and k y than along the {110} 
directions. This suggests that the structural orientation 
of the HFS may be determined by the anisotropy of these 
Fermi surface sheets^ - — 

On the other hand, it may be expected that the strong 
order parameter anisotropy is more important at high 
fields than Fermi surface anisotropy. Since the Kogan 
nonlocal theory does not account for any superconduct- 
ing gap anisotropy, we now consider how this anisotropy 
may influence the VL structure. Both the high field tran- 
sition and the tendency towards a square VL structure 
may be interpreted as caused by the increasing impor- 
tance with field of the predominantly d-wave supercon- 
ducting gap. This is a common prediction made using 
other nonlocal model s 18 ' 19 ' 72 " — and more detailed mi- 
croscopic theories ! 17 ' 22 ' 23 According to the predictions 
of microscopic theory ; 17 ' 22 ' 23 the T of the high field 
square structure is determined by the fourfold vortex 
core anisotropy arising from the d-wave order-parameter. 
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The common prediction of all theories is that the vor- 
tex nearest neighbors of the square VL are expected to 
lie parallel to the nodal direction of the order parame- 
ter JI-i^^iZa-Ts YBasCuaOr at zero field, phase- 
sensitive measurements show that the nodes are not at 
45° to the crystal axes (as expected for a pure d x 2_ y 2 
superconductor), but at ±50° to the b-axisiS as shown 
in Fig [7] (b). The smooth variation of the angle v in the 
HFS phase towards the square value of 90° could reflect a 
field-driven change in the nodal positions of the gap func- 
tion. For the case of YBa2Cu30y, this might indicate a 
change of the admixture which causes the deviation of 
the gap symmetry from pure d x i_ y i. This hypothesis 
could be tested via high-field STM measurements. 

While we may make qualitative comparisons between 
our structural data and the available theoretical work, 
a full understanding of our observations requires inclu- 
sion of the orthorhombic system symmetry, and dis- 
entangling the expected effects due to the anisotropics 
of both the band structure and the nodal gap. This 
is possible in materials where both the band structure 
and gap anisotropies are expected to yield VL structures 
of different orientations, e.g. in Lai.83Sro.i7Cu04+,5J&22 
In YBa2Cu3C>7 however, these two anisotropies are ex- 
pected theoretically to yield similar VL structures! 17 ' 23 A 
more complete understanding of the new structural data 
presented here calls for extensions to the currently avail- 
able theoretical work. In particular, our discussion of the 
Kogan nonlocal theor y 15 ' 16 indicates that instead of re- 
lying on approximate approaches, a better proposition is 
to attempt full first-principles numerical calculations. 



B. Vortex Lattice Form Factor at 2 K 



1. Field- dependence of the form factor 

Within the mixed state, the local field is expressed as 
a sum over its spatial Fourier components at the vari- 
ous different scattering vectors q belonging to the two- 
dimensional reciprocal lattice of the VL. The form factor 
at a wavevector q is the magnitude of the Fourier compo- 
nent F(q), and its value is obtained from the integrated 
intensity, I q , of a VL Bragg reflection as the VL is ro- 
tated through the diffraction condition. I q is related to 
the modulus squared of the form factor, |F((7)| 2 , byZ 6 - 



I q =2 



(4) 



Here, V is the illuminated sample volume, <f> is the inci- 
dent neutron flux density, X n is the neutron wavelength, 7 
is the magnetic moment of the neutron in nuclear magne- 
tons (= 1.91), and $ = h/2e is the flux quantum. Eq. 2] 
shows that the modulus form factor, |-F(<z)|, for a specific 
Bragg spot is obtained from both I q and the magnitude 
g, with all other variables constant. In order to record I q 
experimentally, the sample and cryomagnet are rotated 
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FIG. 8. The angular dependence of the diffracted intensity 
(rocking curve) for the top and bottom Bragg spots of the 
diffraction pattern shown in Fig. [T](a). The dashed lines are 
fits of Lorentzian lineshapes to the data, and the solid lines 
indicate the full width at half maximum (FWHM) . 



together through an angular range that moves a recip- 
rocal lattice vector through the Bragg condition at the 
detector. Fig. [5] shows representative scans of the angu- 
lar dependence of the diffracted intensity (rocking curves) 
for the top and bottom spots with q || a* that form part 
of the diffraction pattern shown in Fig. [T] (a). The quan- 
tity I q is obtained by integrating the area underneath 
the lineshape that is used to fit the rocking curve. When 
appropriate, I q is corrected by the cosine of the angle 
between the scan direction and q (the Lorentz factor) 21 

In Fig. 13 we show our measurements at 2 K of the VL 
form factor across the entire field range up to 10.8 T. At 
each field, we distinguish between form factor values ob- 
tained from different types of Bragg spots, and each value 
represents a mean taken over equivalent spots. Close to 
the first-order VL structure transitions, the form factor 
values are noticeably lower than those expected from an 
extrapolation of 'pure phase' data recorded away from 
these field regions. This is owing to the strong field- 
dependence of the occupation fraction of each of the two 
relevant structures close to the transitions, and is an ef- 
fect not accounted for by Eq. 2J For the low field tran- 
sition, one sees reduced form factor values only at fields 
of 1.75 T and 2.5 T. For fields within this range, we were 
unable to measure reliably the integrated intensities for 
the Bragg spots due to a disordering of the VL. For the 
high field transition, reduced form factor values are seen 
at 6.5 T and 7 T. In this case, the two structures co- 
exist cleanly and it is possible to account for all of the 
diffracted intensity. However, when modeling our form 
factor data, we only include data obtained at fields where 
there is just one VL structure type in the sample. 

The data shown in Fig. |9] indicate that with increas- 
ing field, an increasingly significant anisotropy emerges 
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FIG. 9. (Color online) The field-dependence at 2 K of the 
observed form factors f(q). Across the entire field range, 
we distinguish between the different types of Bragg spot ac- 
cording to the different magnitudes of the q- vectors for each 
type. Each datapoint represents the average form factor for 
a certain spot type. 



FIG. 10. (Color online) The field-dependence at 2 K of the 
ratio between on-axis and off-axis form factors. For the LFS 
and HFS phases, the on-axis spots lie along a*, whilst for 
the IFS phase, the on-axis spots lie along b*. The black 
dashed line indicates a ratio of unity, which is expected in the 
anisotropic London limit. 



between the form factors of the two types of Bragg spot. 
This is explicitly shown in Fig. 1101 where for each VL 
structure type we present the field-dependence of the 
form factor ratio |.F(q on-axis) | / |-F(q off-axis) |. We see 
that this denned ratio varies monotonically over the much 
of the field range, even though the primitive cell of the 
IFS phase is orthogonal to that of the LFS and HFS 
phases. For the LFS phase, the data in Fig. [10] show 
that all the spots have the same form factor within the 
experimental uncertainty, thus giving a form factor ratio 
close to unity. This is as expected from London theory, 
even for the case of an anisotropic VL. In moving from 
the isotropic case to the anisotropic case appropriate for 
YBa 2 Cu307, F (q) changes according to 



F(q) 



(B) 



(B) 



1 + q 2 X 2 1 



(5) 



Here, (B) is the average internal induction (which 
within experimental error is always = due to the 

plate-like crystals and large n) and A^ is the penetration 
depth along the directions The terms q x and q y rep- 
resent the x and y components of the q- vector as viewed 
in the diffraction patterns shown in this paper; thus q x 
is parallel to b*, and q v parallel to a*. In this case, it 
may be shown^ that the a-b anisotropy in the q- vectors 
in this distorted hexagonal structure exactly cancels the 
penetration depth anisotropy, giving the observed result. 
Another prediction of London theory is that the form 
factor is independent of field for B B c \ (where the 
1 in the denominators of Eq. [5] may be neglected). The 
fall-off shown in Fig. [5] will be discussed in the follow- 
ing section. For the IFS phase, the observed form factor 



anisotropy unusually has \F(q || b*)| < |F(q off-axis)| 
giving a smaller value of for the Bragg spot with 

the shorter q-vector. 



2. Modeling the field- dependence of the VL form factor 

The simplest analytic models that are most commonly 
used to interpret form factor data are those based on ap- 
plying corrections to local electrodynamics suggested by 
Ginzburg-Landau (GL) theoryiSr— Although GL theory 
is only strictly valid close to T c , the results of more ex- 
act numerical calculation o 22 ' 82 ' 83 show that such simple 
models^r— provide a reasonable description of the mixed 
state at lower T . In the isotropic case, the model is 



2t2\ 



F{q) = 



(-B)exp (— eg £' 
1 + g 2 A 2 



(6) 



which can be considered as a refined London model with 
a gaussian cut-off term that accounts for the finite size of 
vortex cores. Here, £ is the in-plane GL coherence length 
and A is the in-plane penetration depth. The constant c 
is a parameter expected to have a value between 1 /4 and 
2i£i A quantitative comparison between the predictions 
of the above London model and those of numerical calcu- 
lations carried out within the Eilenberger formalis m 22 i 84 
suggest that over our field range, and at low T, an ap- 
propriate value for c is 0.44^ We note that this value 
lies close to that of |, which has been used in successful 
form factor analyses for other superconductorsi 52 i 54 

Other authors have used the Clem model^ which is 
obtained from a variational solution to the GL equations. 
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TABLE I. A summary of the returned parameters obtained after fitting the 2 K form factor data with both the London and 
Clem models (see text for details of the fitting). The values of \ 2 are those obtained after each step of the analysis procedure. 



Model 


\ a (nm) 


\b (nm) 


£a (nm) 


£t> (nm) 


X 2 (Step 1) 


X 2 (Step 2) 


London 


138(2) 


107(1) 


3.04(4) 


3.54(11) 


1.46 


1.92 


Clem 


132(5) 


103(1) 


2.72(13) 


2.49(36) 


4.40 


13.27 



The model is: 



F(q) = (B) 



9K1 (g) 

1 + q 2 X 2 ' 



fJ 



V2£ 



GL (q 2 + \- 



>V2 



(7) 



where K\ is a modified Bcsscl function of the second kind. 
This model appears advantageous in that there is no ill- 
defined c factor within the core-correction term gKi(g) 
as there is in Eq. [51 but it relies on the GL equations 
which are only valid close to T c . As we will show, this 
approach gives a much poorer fit to our data than the 
refined London model described by Eq. [5] 

Both of these models can be extended in the same 
way in order to account for the biaxial anisotropy of 
YBa 2 Cu307. For example, the refined London model 
(Eq. [6J can be re- written in the convenient form, 



^(q) 



(B)exp (-0.44 + 



q 2 A 2 



(8) 



Here £j represent the GL coherence lengths along the di- 
rections i. According to the extended models, in general 
the form factor is expected to depend on the four param- 
eters Aj and However, the number of free parameters 
can be reduced by fitting the model independently for 
certain types of spots. For spots with q |j a*, terms de- 
pendent on q x vanish, as do those dependent on q y for 
spots with q || b*. Taking advantage of this, we extract 
values for Ai and £j from our form factor data using the 
following analysis procedure: 

Step 1. We fit the form factor data for Bragg spots with 
q || a* in order to obtain A& and £ Q . 

Step 2. The values of A& and £ a found in Step 1 are held 
constant, and we fit the data for the spots with an 
off-axis q- vector in order to obtain £ a (0) and Ab (0). 

For both steps, the values of q x and q v at each field, and 
their associated errors, are those determined experimen- 
tally. We also see that the fit in each step is sensitive to 
just two fitting parameters. 

In Fig. QT] we show the fits of the form factor data to 
the London model. The first step of the analysis pro- 
cedure is shown in the inset, and yields parameters of 
A5 = 107(1) nm and £ a = 3.04(4) nm. The reasonable 
fit quality is indicated by x 2 — 1.46. In the main panel 
we show the results of the second fit, with the fitted pa- 
rameters A a = 138(2) nm and = 3.54(11) nm, and 
X 2 = 1.96. As can be seen in Table|H using the same anal- 
ysis procedure with the Clem model provides a poorer de- 
scription of the data than the London model. Attempts 
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FIG. 11. (Color online) A semi-logarithmic graph showing the 
field-dependence of the VL form factor for Bragg spots with 
a q-vector lying off-axis. The fit of the data to the London 
model is also shown. The inset shows the model fit of the 
form factor data with q || a*. 



were made to fit the data using a version of the Clem 
model that includes minor corrections considered appro- 
priate for high-T c superconductors^, However, these did 
not generate an improvement of the fit. 

We now discuss the results obtained on fitting the 
form factor data to the London model (Eqn [5]). The 
values returned for A a and Af, are in reasonable agree- 
ment with those determined by previous SANS studies 
on optimally- and over-doped YBa2Cu307_,$r£ though 
overall our values are a little shorter. We find better 
agreement between our values and those found in a recent 
/iSR study on a clean sample of close to optimally-doped 
YBa 2 Cu306.92^ This indicates the state of our sample 
to be of higher quality and more completely oxygenated 
compared to the sample studied previously by SANS<2& 
On the other hand, the values returned for £ a and £b are 
likely to be overestimates. Using these values we can esti- 
mate the upper critical field to be of order ~ 30 T, which 
is clearly too low for YBa2Cu307. One possibility for 
the over-estimation of £j is that VL disorder may reduce 
the measured intensity in the form of a low T 'static' 
Debye- Waller (DW) factor. Such a factor would enter 
the numerator of Eq. [8] as exp (— g 2 (u 2 )/4), where (u 2 ) 
is the root-mean-square displacement of a vortex along 
the direction of q. By inspection it is clear that this term 
has precisely the same effect on the VL form factor as the 
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Gaussian cut-off term that accounts for the finite vortex 
core size, and that both of these terms become increas- 
ingly important at larger q. As a consequence of this, we 
are unable to disentangle the effects of £ or (u 2 ). How- 
ever, bearing in mind that in high-T c materials both the 
value of £ is short, and the vortices are 'soft' (small val- 
ues of the shear modulus and the tilt modulus at small 
distances 87 ) it is unsurprising that disorder in the form 
of a DW factor may become important. This is partic- 
ularly the case at higher fields where the VL spacing is 
just a few times the value of £j . 

An alternative possibility is that any of £j or A^ exhibit 
a significant field-dependence over our field range. Previ- 
ous high field fiSR studies on YP^CusOe.gs report both 
a strong increase with field of the extracted in-plane pen- 
etration dept h 27 ' 28 attributed to non-local effects ) 20 ' 88 
and also a strong decrease with field of the vortex core 
size £ attributed to a field-induced quenching chain su- 
perconductivity. — It should be noted that in the /zSR 
analysis 28 an undistorted hexagonal VL was assumed, 
which we have shown is not the case. Also, the Clem 
model was used in a T range where its validity is ques- 
tionable. Within the framework of our form factor model, 
unconventional field-dependences of \ and £j may offset 
one another to give alternative fits to the data. However, 
we hesitate to fit for field-dependent length-scales as this 
amounts to inserting unknown free parameters. 

In spite of our conclusion that the values for £j are 
too large, the returned values for A a = 138(2) nm and 
As = 107(1) nm are likely to represent the intrinsic bulk 
values. Within the local London model described by 
Eq. [51 the values of Aj are determined by the form factor 
magnitude at lower fields, and are hence rather insensi- 
tive to any exponential terms in the numerator. From 
our analysis we find that — \ a /\ — 1-29(2), a value 
which is in good agreement with that obtained directly 
from the measure of the distortion of the hexagonal struc- 
ture shown in Fig. [3J This strongly indicates that local 
theory is at least able to describe the anisotropy of the 
low field VL. In addition, the evolution of the form factor 
anisotropy with field [Fig. [10] could arise, at least in part, 
if the value of j\ does not fall significantly with increas- 
ing field. In this case, the form factor anisotropy emerges 
as a consequence of the changes in VL structure and the 
associated changes in q, rather than due to anisotropic 
field-induced variations in the fundamental length-scales. 

We emphasize that conclusions drawn from form factor 
analyses such as ours, and also those of /iSR data, need 
to be considered with great care, and provide only an 
approximate picture of the mixed state. A more thorough 
understanding of our form factor data might be obtained 
by harnessing the power of a microscopic theor y 21 ! 22 ' 84 
where additional effects such as those due to non-locality 
and order-parameter anisotropy are included by default. 



C. Measurements at Higher Temperatures 

1. VL Structure at Higher Temperature 

Our T-dependent measurements were carried out by 
warming from, or cooling into, an oscillation field-cooled 
state. While the T was changing, the field was oscillated 
using the OFC procedure. Upon reaching the desired T, 
the field was held stationary at the intended value be- 
fore starting the measurement. In Fig. [T^] we present 
T-dependent measurements of the VL structure open- 
ing angles defined in Fig. Q] at various magnetic fields. 
We checked for hysteretic behavior by performing both 
warming and cooling T scans at 1.0 T [Fig. [T^] (a)] and 
5.0 T [Fig. (b)]. At these two fields, there is no dis- 
cernible difference in the precise VL structure between 
warming and cooling, thus showing that hysteresis effects 
on the VL structure can be neglected in our sample. 

Fig. [12] (a) reveals that for fields within the LFS phase, 
the VL structure remains essentially T-independent. 
Only when T — > T c is there any measurable variation in 
the opening angle </>, though this variation remains small. 
Data obtained in the IFS phase [Fig. 1121(b)] indicate that 
this weak T-dependence of the VL structure persists un- 
til fields of up to 5.0 T. However, more pronounced In- 
dependences emerge at higher fields. At 6.9 T, a field 
where the IFS and HFS phases co-exist at 2 K, there is 
clear T- induced reduction in the values of both p and v as 
T — » T c . At yet higher fields in the HFS phase, the data 
at 8.0 T show the high T reduction of v to become in- 
creasingly pronounced. Remarkably, the data in Fig. Q21 
reveal that for no field does adjusting the T cause a tran- 
sition between different VL structures. Even at 6.9 T, 
the data show that the low T co-existence of the IFS and 
HFS phases persists over the entire T range. All of these 
observations suggest that the phase boundaries between 
different structure types are unusually T-independent. 

It also remains to be explained why any variation in 
the VL structure occurs only at high T, and at relatively 
high fields. The former observation suggests that, in spite 
of utilizing the OFC procedure, the VL structure is still 
'frozen-in' and pinned at a relatively high T m of ~50 K. 
If so, any variation of the VL structure with T would 
only occur in the reversible region, when the vortices are 
depinned and able to adopt the equilibrium structure. 
For fields below ~ 5.0 T, there is little variation of the 
structure at higher T, suggesting that there will not be 
a significant difference between the precise VL structure 
frozen-in at Tj rr , and that expected for low T . However, 
at higher fields a larger difference emerges between the 
VL structure frozen-in, and that which may be expected 
at lower T by extrapolating the high T behavior. 

The high T variation of the VL structure can be un- 
derstood within the framework of non-local theory 15 ' 17 ' 19 
where the non-locality can arise either from Fermi surface 
anisotropy or gap structure. A common prediction is for 
non-local effects to become suppressed as T — > T c . This is 
most clear to see at fields within the HFS phase, where on 
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FIG. 12. (Color online) The temperature-dependence of the VL structure opening angles (a) <j> and v for the LFS and HFS 
phases respectively, and (b) p for the IFS phase. The indicated angles correspond to those in Fig. [T] In both panels, the lines 
passing through for the data at each field are guides to the eye. Data are shown only for temperatures where the spot position 
on the detector could be determined accurately. Error bars not visible can be considered of order the size of the datapoint. 



the approach to T c we observe the values of v to approach 
those values of <j> observed in the LFS phase where we 
believe that anisotropy due to local effects is dominant. 
Non-local effects are also suppressed by T within the IFS 
phase, and the high-T variation in p shows that the VL 
tends to become slightly more anisotropic (larger value of 
rj) close to T c . In this case, the suppression of non-locality 
allows the VL to adopt a hexagonal coordination that is 
more distorted. A common interpretation of these data 
is that the high T or low field anisotropy factor (~1.29) 
just arises due to the band mass anisotropy, and can be 
understood with local anisotropic London theory. 



We now consider the remarkable T- independence of the 
VL structure transition fields. It is expected that non- 
local effects fall off with increasing T but not to zeroi 15 ' 19 
Clearly at T—T C (H), these effects remain strong enough 
to control the orientation of the VL primitive cell, and 
it appears that the crossover fields remain essentially T- 
independent below this. However, in the IFS and HFS 
phases the shape of the VL evolves away from square 
structures with increasing T. This tendency is also con- 
sistent with that expected if the VL structures are stabi- 
lized by the (i-wave gap anisotropy 22 ' 75 With increasing 
T, the influence of the order parameter anisotropy on the 
VL structure is expected to be suppressed, as quasiparti- 
cles that previously occupied states near the nodal posi- 
tions adopt a greater angular spread in reciprocal space. 
However, if the d-wave gap anisotropy is responsible for 
driving the VL structure transitions, and controlling the 
primitive cell orientation, the influence of this anisotropy 
must persist to surprisingly high T. 



2. VL Form Factor at Higher Temperatures 

It seems that the precise VL structure is frozen-in 
on cooling below Tj rr ; this provides the opportunity to 
record the T-dependence of the form factor at low fields 
without requiring an excessive amount of neutron beam- 
time. Typically, carrying out these measurements is a 
lengthy process; at each T, full rocking curves such as 
those shown in Fig. [8] should be obtained for each Bragg 
spot. However, in our case where the precise VL struc- 
ture is weakly T dependent, the same information can 
be obtained by counting solely at the peak of the rock- 
ing curve. Knowing this peak intensity, along with the 
FWHM of the rocking curve and the q-vector allows a 
calculation of the form factor. This approach assumes 
that the FWHM of the rocking curve remains T indepen- 
dent. This assumption was checked within each T scan. 
While at most temperatures, measurements were carried 
out solely at the Bragg angle, at certain temperatures 
we carried out full rocking curve measurements. These 
measurements provide both a check of the assumed T 
independence of the rocking-curve FWHM, and also full 
measures of the form factor to be compared with those 
obtained from the fixed angle measurements. 

In Fig. [13] we present the T dependence of the form 
factor in applied fields of 0.2 T (Fig. [TJ (a)) and 5.0 T 
(Fig. [T2] (b)). At each field, and for each type of Bragg 
spot, we show both the values obtained by measuring at 
fixed rotation angle and those obtained with the usual 
full rocking curves. The full rocking curve measurements 
allow us to plot the T dependence of the FWHM of the 
rocking curve for each type of Bragg spot. Both the insets 
to Figs. [13] (a) and (b) show that there is no clear sys- 
tematic T variation. Therefore, we use a mean value for 
the FWHM of each type of Bragg spot across the entire 
T range. As seen in Figs.[T3](a) and (b), the agreement 
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FIG. 13. (Color online) The temperature dependence of the form factor |T(g)| for the two types of Bragg spot in (a) the LFS at 
0.2 T, and (b) the IFS at 5.0 T. For each scan, we performed an OFC procedure down to base temperature, and then measured 
on heating the sample. In (a), the solid and dashed lines show fits to the data made using both a minimal and local d-wave 
model. In (b), solid and dashed lines show fits to the data after extending the local d-wave model used in (a) with a non-local 
correction. Both models are described in the text. The inset of each figure shows the temperature-dependence of the rocking 
curve FWHM as obtained from the rocking curve measurements. The dashed and dotted lines, respectively, indicate the mean 
value of the FWHM used to normalize the form factor values of the fixed angle data (see text for details). 



between the two techniques of obtaining the form factor 
is good, thus providing confidence in our measurements 
obtained at fixed rotation angle. We also note that any 
T-induced variation in the precise VL structure is small, 
and does not extend beyond the in-plane resolution of 
the instrument. Therefore, at these fields, any change in 
the VL structure will have a negligible influence on the 
values of intensity obtained at fixed angle. 

In Fig. Q2] (a) we see that at 0.2 T the form factors 
for the two types of Bragg spot exhibit a similar In- 
dependence. Interpreting these data requires a knowl- 
edge of the T-dependences for both £j (T) and Xi (T), 
each of which require a detailed knowledge of the band 
structure. While ARPES studie a 89 ' 90 and various calcu- 
lations have established the main band structure reason- 
ably well ) 70 i 71 computing the T dependence of the super- 
fluid density in materials like YBa2Cu30y is not trivial 
since this requires a precise knowledge of the momentum- 
resolved quasiparticle spectrum. As this information is 
not readily available, we instead make use of a 'min- 
imal model' for the Fermi surface sheets arising from 
the Cu02 planes, and represent them by a single two 
dimensional quasi-cylindrical sheet, which has some a- 
b anisotrop y 70 i 71 ' 89 i 90 Whilst such an approach neglects 
a possible contribution to the superfiuid density due to 
chain superconductivity, it is expected that the dominant 
contribution will arise due to the large hole-like sheets 
centered at the Brillouin zone corner. 

According to Eq.[51 the form factor F (q, T) for a Bragg 
spot will be temperature-dependent according to: 



p i rp\ (j?)exp(-0.44(q^ 2 (T) + q^(r))) 



In the local approximation, the T-dependence of the pen- 
etration depth can be calculated using the following equa- 
tion suitable for a two-dimensional and cylindrical Fermi 
surface i21 



1 



K,b (T) 



Ps,a/b (T) 



4ttT 



cosh 



2k B T 



d(f>de. (10) 



Here, 4> is the azimuthal angle around the cylindrical 
Fermi surface and \/e 2 + A 2 (T, <fi) defines the excitation 
energy spectrum. The gap function A k (T, <j)) is assumed 
separable into momentum- and T-dependent factors such 
that A k (T, <j>) = A k {(p) A (T) where A k (0) describes 
the angular-dependent part of the gap function, and 
Ao (T) the T-dependent part. The angular dependence 
for the gap function can be expressed as A& (<f>) — 1 for an 
s-wave gap, or A k (<fi) — cos(20) for a c£-wave gap. The 
T-dependence of the gap is calculated using the following 

expression, 92 A (T) = A (0) tanh (l.78y/T c /T - l) , 

where Ao (0) is the maximum magnitude of the gap at 
zero T. This gives a good approximation to the BCS 
weak-coupling T-dependence, though A does not need 
to have the BCS value. We also use this expression for 
the gap to calculate the T-dependence of ^ a ,b (T), 



ia/b (T) = U/b (0) tanbVl.78 (T c /T - 1) 



(11) 



q 2 A 2 (T) + q 2 A 2 (T) 



By using the values of & and \ obtained from the 2 K 
form factor analysis as the zero temperature values, the 
only fitting parameter for a single gap model is the zero 
temperature gap magnitude Aq (0). 
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By making use of a single d-wave gap function, 
Fig. Q2] (a) shows that the minimal model is able to pro- 
vide a good description of the data at 0.2 T. The returned 
values for Ao (0) are 25(2) meV for the Bragg spot with 
q || a*, and 26(2) meV for the off-axis Bragg spots. These 
values are in reasonable agreement with those found in 
/iSR experiments on twin-free YBa2Cu307_a. 51 

On the other hand, and again drawing comparison be- 
tween our results and those in the aforementioned /iSR 
study* 5 - 1 - we do not observe a signature of an inflection 
point in \F (q,T)\ at low T. This inflection point, which 
is observed clearly in directly comparable low field mea- 
surements of the superfluid density is attributed to a con- 
tribution to the superfluid density of a small s-wave gap 
that is quenched with increasing field. 51 An absence of an 
inflection point in microwave studies has previously been 
interpreted as possibly attributable to chain disorder, 46 
which can not be the case in our fully oxygenated sample 
of YBa2Cu30y. We also stress that even though the VL 
is pinned at low T, an increase of the superfluid density 
would be clearly visible in our measurements of the form 
factor at 0.2 T. 

In contrast to the data obtained at 0.2 T, in Fig.[T3"lfb) 
we see very different behavior at 5.0 T. For each spot 
type, the low T-dependence of the form factor is much 
weaker, and the minimal model is less able to provide a 
description of the data. To obtain a better agreement 
we would need to invoke either a vortex core contraction 
with increasing T, or a T-dependence to the superfluid 
density that is more reminiscent of fully-gapped behavior. 
Both of these scenarios are considered as unlikely. 

An alternative route to resolving this issue is to realize 
that a weak low-T dependence of the form factor is in bet- 
ter agreement with that expected due to a strong nonlo- 
cal response of a <i-wave superconductor under an applied 
field32&22 Amin et a/. 20 ' 88 have studied the effect of a 
field-induced nonlocal response on the effective penetra- 
tion depth in YBa2Cu306.95, as observed by //SR. 28 With 
increasing field, they find that the usual linear low-T de- 
pendence of the superfluid density crosses over to a T 
dependence below a temperature T* = Ao(£o/^) oc y/~H, 
where d is the VL spacing. These results suggest that 
the weak low T-dependence of \F (q, T)| we observe may 
also be attributable to strong nonlocal effects. 

Consistent with the theory of Amin et a/. 20 ' 88 , a simple 
expression that captures the expected behavior of the 
superfluid density due to d-wave non-local effects is 

= 1 - (1 - ?.,„/» CO) (^jf-) ( T , ^ T . f ) (12) 

where p s ,a/b IS the superfluid density as calculated in the 
local limit (Eg. ITD|) . and n s a /i, is the superfluid density 
when accounting for the non-local effect. By combining 
Eqs. I^and [T21 it is possible to calculate the temperature 
dependence of the form factor with a non-local correction 
to the penetration depth. 



A fit to the data was most easily achieved by using 
pre-determined values for A (0) and T*, and leaving all 
of & and Xi as free parameters. A (0) was fixed to be 
25.5 meV, the mean value found for the 0.2 T data. The 
value of T* was held at the temperature where linear 
extrapolations of the low and high temperature portions 
of the dataset intersect. At 5.0 T, and for both spot 
types, this intersection temperature is 51.4(5.0) K. Fur- 
ther analysis showed that the quality of the fit is actually 
quite insensitive to T* , and as such its precise value does 
not significantly affect the outcome. 

Similar to the 2 K form factor analysis, the data were 
fitted using a two step analysis sequence: 

Step 1. We fitted the data for the Bragg spot with q || b* 
in order to obtain values for £j (0) and A a (0). 

Step 2. The values of £& (0) and A a (0) found in Step 1 
are held constant, and we then fit the off-axis data 
in order to obtain the values of £ a (0) and (0). 

Fig. [13] (b) shows the good fits of the nonlocal model to 
the data at 5.0 T with the values obtained from the fitting 
procedure being A a (0) = 170(4) nm, A 6 (0) = 121(3) nm, 
£ Q (0) = 1.82(18) nm and (0) = 1.55(15) nm. However, 
a good fit is only obtained using the nonlocal model when 
all of £j (0) and Xi (0) are left as free parameters. Other 
approaches were tried: using the local model (Eq. [UJ) with 
field-dependent length-scales £j (0) and Xi (0); or using 
the non-local model, but keeping the values of £j (0) and 
Aj (0) fixed to those obtained from the 2 K form factor- 
analysis. Both of these gave significantly poorer fits. 

By comparing the values of £i(0) and A^(0) obtained 
from fit shown in Fig. [T3] (b) , to the values obtained from 
the field-dependent form factor analysis at 2 K, we might 
conclude that A^(0) increases with field, while £i(0) de- 
creases with field. The signs of these field-induced vari- 
ations are consistent with those reported by fj.SH stud- 
ies on YBa2Cu3 06.95' 27l2S The possible physical reasons 
for these field-dependences were introduced previously 
in our discussion of the field-dependence of the form fac- 
tor at 2 K. As was also mentioned in that discussion, 
before making firm statements on any field-dependence 
for either of £i(0) or Aj(0), the possible role that VL 
disorder may play on our measurements should be con- 
sidered. Since our measurements of the rocking curve 
FWHM shown in the inset of Fig. Q]|] (b) do not exhibit 
a large T-dependence, we may conclude that a term in 
the form of the static DW factor does not dominate the 
low T-dependence of the VL form factor for wavevectors 
associated with fields of up to 5.0 T. Nevertheless, more 
detailed measurements may reveal a slight T-dependence 
of the rocking curve FWHM, and thus provide alterna- 
tive descriptions of the data than those solely based on 
field-dependent length-scales. At 0.2 T, ignoring disorder 
is more easily justifiable, because the DW factor is closer 
to unity at the shorter wavevectors associated with low 
fields. In principle, the role of VL disorder can be deter- 
mined from direct measurements of higher-order spots of 
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FIG. 14. (Color online) The temperature-dependence of the 
form factor ratio taken between the different types of Bragg 
spots at at various fields. For the fields of 0.2 T and 8.0 T, 
taken in the LFS and HFS phases respectively, the form factor 
ratio plotted at each temperature is \F(q || a*)| / |-F(q |fa*)|. 
For the 5.0 T data taken in the IFS phase, the ratio is 
|F(q||b*)|/|F(qrb*)|. 



the VL because these spots are more sensitive to disor- 
der effects. Unfortunately, due to the long value of A in 
YBa2Cu307, investigating higher-order spots is experi- 
mentally challenging. Indeed, within the IFS phase, we 
were unable to observe any signal due to higher-order 
spots of the hexagonal VL. 

Despite the ambiguity caused by possible VL disorder, 
a term in the form of the DW factor cannot explain the 
emergence of the low T form factor anisotropy between 
the different types of Bragg spot. Instead, this anisotropy 
is likely to arise as a consequence of field-induced non- 
local effects on the VL structure. This is further inferred 
by an examination of the T-dependence of the form fac- 
tor ratio. Data for selected fields are presented in Fig.fTTl 
where, as for Fig. [TU1 at each T we plot the form factor 
ratio between the on-axis and off-axis spots. Within the 
LFS phase, the data at 0.2 T show the ratio to remain at 
unity over the entire T range, which is as expected for a 
VL close to the local London regime. At 5.0 T the form 
factor ratio, which is less than unity at low T, smoothly 
tends towards unity on the approach to T c . These re- 
sults show that increasing T suppresses the low T form 
factor anisotropy that is induced by nonlocal effects on 
the VL. The persistence of this behavior to fields within 
the HFS phase is also tentatively supported by our re- 
sults obtained at 8.0 T. This indicates that it may be 
possible to move between a square-like phase with four 
strong spots and two weaker spots (a VL like that shown 
in Fig.[T](c)) towards a hexagonal-like structures with six 
spots of similar intensity (like that shown in Fig. [TJ (a)) 
just by increasing T. We expect future SANS measure- 
ments to shed light on the validity of this suggestion. 



FIG. 15. (Color online) The (^ H,T) VL structure phase di- 
agram for fields applied parallel to the c-axis. Circle data 
points were obtained from measurements of the VL form fac- 
tor, and correspond to the point where each relevant VL struc- 
ture is measured to occupy 50% of the sample volume. Di- 
amond datapoints correspond to estimates of the transition 
points as determined from structural measurements. The dot- 
ted and dashed lines are guides for the eye phase boundary 
lines determined from the form factor data only. The melting 
line is deduced from data presented in Ref. (9J. 



3. Phase Diagram for fioH || c 

In this final section, we summarize the main results 
of our study in the form of the fioH | c phase diagram 
shown in Fig 1151 Datapoints were determined from mea- 
surements of either the VL structure or the VL form fac- 
tor. For the LFS - IFS transition, our best estimate of the 
transition field at 2 K using form factor data is 2.3(2) T, a 
field slightly higher than that determined from the struc- 
tural data alone. For the IFS - HFS transition at 2 K, the 
transition field determined by both structural and form 
factor measurements are in good agreement, both being 
at 6.7(2) T. This phase diagram shows more clearly how 
the phase boundary lines between different VL structure 
phases remain essentially constant in field across the en- 
tire T range, a result consistent with that inferred from 
measurements of the T-dependence of the VL structure. 
While we can not rule out any curvature in either phase 
boundary line closer to T c , the SANS measurements show 
no evidence for a strong variation. 

We also draw comparison between the phase diagram 
shown in Fig [151 and that obtained from SANS stud- 
ies on a lightly twinned sample of YBa2Cu307. 39 The at 
best weakly T dependent phase boundaries observed in 
the twin-free sample are in contrast to the single phase 
boundary line observed in the twinned sample, which 
clearly curves upwards in field at higher T£2- Further- 
more, as the structure phase boundary for the twinned 
sample is continuous, we suggest that its field and T 
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phase diagram is determined by a complex balance be- 
tween the intrinsic T of the high field VL, and the ex- 
trinsic T associated with the {110} pinning potential. 



IV. SUMMARY AND CONCLUSIONS 

In this study we have reported a small-angle neutron 
scattering (SANS) study of the VL in a high quality, and 
twin-free, sample of YBa2Cu3C>7. The absence of vor- 
tex pinning to twin planes allows the intrinsic VL struc- 
ture to be imaged, and provides the possibility to make 
more direct comparisons between theoretical prediction 
and field- and temperature (T)-dependent observations 
than is otherwise possible in twinned samples. 

Our measurements of the VL structure at 2 K show 
that a field-driven sequence of first-order structure tran- 
sitions exist in twin-free YBa2Cu30y. The driving mech- 
anisms behind these transitions can be broadly described 
as caused by the increasing importance with field of non- 
local effects on the intervortex interaction. Theoretically, 
non-local effects are catered for by the development of 
higher-order correction terms to local theories. These 
terms play the role of coupling the VL to anisotropies of 
the host material. One source of anisotropy expected to 
influence the VL properties is that due to the band struc- 
tureji^~— with another source arising due to the predom- 
inantly d x 2„ a 2-wave order parameter anisotropy . 18 i 19 

Within the theoretical literature, we find that no model 
is able to provide a full description of the field-dependent 
sequence of VL structure transitions that we observe. 
The simplest theory that seems most appropriate for 
making predictions of the field-dependent VL structure in 
twofold symmetric systems like YBa2Cri3 07 is provided 
by Kogan et aZj 15 i 16 This model accounts for nonlocal 
effects that arise due to the band structure anisotropy, 
and not superconducting gap anisotropy. Our close ex- 
amination of the Kogan theory has revealed some hith- 
erto unreported properties which show that in its current 
form, the model framework is over-simplified and unable 
to explain field-dependent observations in real twofold 
systems. Nevertheless, in spite of the noted shortcom- 
ings of the model, there are indications within our data to 
suggest that nonlocal effects that originate from the band 
structure anisotrop y 15 i 16 may be the cause of our low field 
transition between orthogonal hexagonal VL structures. 
Similar arguments can be made to suggest that the high 
field transition between hexagonal and rhombic VL struc- 
tures has the same origin. However, this second transi- 
tion can also be understood as caused by the increasing 
influence of the predominantly ci-wave gap anisotropy, 
whether this tendency is captured by a non-local the- 
or y 18 ' 19 or by detailed microscopic calculations! 17 ' 22 ! 23 
Based on our structural data alone, it is difficult to de- 
cide which of these two sources of anisotropy causes the 
high field square-like structure to appear. The available 
predictions lead us to expect a high field square-like VL 
structure with an orientation consistent with that which 



we observe, regardless of whether band structure or pair- 
ing anisotropy is dominant. However, measurements at 
higher T show each of the structure phase boundary lines 
to exhibit the same weak T-dependence, suggesting that 
a common mechanism may lie behind both transitions. 

In the lowest field VL structure phase, local anisotropic 
London theory^— appears to provide an adequate de- 
scription for the VL distortion and form factor. While 
such a theory is unable to explain the stabilization of a 
preferred orientation (which must be due to weak non- 
local interactions), it does provide a framework within 
which we can obtain a measure of the in-plane anisotropy 
parameter 7a = A a /Af,, either directly from the measure 
of the distortion of the low field hexagonal structure, or 
from fitting parameters obtained in an analysis of the 
field-dependence of the form factor at 2 K. 

An important finding is that a regime close to the lo- 
cal limit ends abruptly on moving into the higher field 
structure phases. While the strongest evidence for high 
field non-local effects is provided by our measurements 
of the VL structure, other evidence is found on study- 
ing the field- and T-dependence of the VL form factor. 
In particular, the emergence of a form factor anisotropy 
that becomes increasingly significant at larger fields, and 
which is quenched with increasing T, points towards 
the decisive role that non-local effects play on the ob- 
served VL structure. Furthermore, within the interme- 
diate field phase at 5.0 T, we observe an unusually weak 
low T-dependence of the VL form factor for both types 
of Bragg spot. An analysis using a simple model that in- 
cludes a non-local correction suggests that this unusual 
behavior may arise due to non-local effects induced by 
the anisotropy of the d-wave order parameteri 20 i 88 While 
such a d-wave non-local theor y 20 ! 88 may explain the ob- 
served weak low T-dependence of the VL form factor at 
5.0 T, there is poor agreement between the VL structures 
that we observe, and those predicted within the same 
theoretical frameworkji 8 . - — This is likely because these 
models do not consider a vortex core anisotropy that 
will arise due to anisotropies in either the band struc- 
ture or the gap function, and which other theoretical 
studies show can be responsible for field-driven transi- 
tions between different VL structures! 17 i 22 i 23 Overall it 
seems that the precise field- and T-dependent VL prop- 
erties are sensitive to a delicate balance between the field- 
dependent anisotropies in the system. The high quality 
data presented in this study should serve to inspire the- 
oretical interest in the precise properties of the VL in 
High-T c YBa2Cu30y. Further progress in understanding 
these results will also benefit from new SANS studies at 
yet higher magnetic fields. 
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